Constraints on the Holographic Dark Energy Model from Type la Supernovae, 
WMAP7, Baryon Acoustic Oscillation and Redshift-Space Distortion 
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In this paper, we use the joint measurement of geometry and growth rate from matter density 
perturbations to constrain the holographic dark energy model. The geometry measurement includes 
type la supernovae (SN la) Union2.1, full information of cosmic microwave background (CMB) 
from WMAP-7yr and baryon acoustic oscillation (BAO). For the growth rate of matter density 
perturbations, the results f{z)as{z) measured from the redshift-space distortion (RSD) in the galaxy 
power spectrum are employed. Via the Markov Chain Monte Carlo method, we try to constrain 
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the model parameters space. The jointed constraint shows that c 

as = 0.763lo o4g5lo og28io 108 with l,2,3cr regions. After marginalizing the other irrelevant model 
parameters, we show the evolution of the equation of state of HDE with respect to the redshift z. 
Though the current cosmic data points favor a phantom like HDE Universe for the mean values of 
the model parameters in the future, it can behave like quintessence in 3cr regions. 
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I. INTRODUCTION 

The holographic principle says that the number of de- 
grees of freedom in a bounded system should be finite 
and has relations with the area of its boundary By 
applying the so-called holographic principle to cosmol- 
ogy, one derives a relation between vacuum density and 
a cosmological scale pA = 3c'^M^iL~'^ [IHll: where c is a 
numerical constant and Mpi is the reduced Planck Mass 
Af-2 ^ gjrG. The obtained vacuum energy, dubbed as 
holographic dark energy (HDE), can push our Universe 
into an accelerated expansion phase at late time [1, Q . 
By taking different cosmological scale, for example the 
Hubble horizon [H 0, Q , the event horizon or the parti- 
cle horizon Q as discussed by and the Ricci scalar 
0, one has different HDE model. Based on the idea 
that gravity as an entropic force Q, a similar DE den- 
sity was given in Q where a linear combination of 
and H was also presented, see also [Hi- Further- 
more generahzed HDE models pR = Sc^ M^iRf {H'^ / R) 
and ph = 3c^MpiH'^g{R/H'^) were also presented in Ref. 
[T^ . In this paper, we consider the typical HDE model 
where the future event horizon 



Reh{a) = a 



dt 



da 



(1) 



is taken as a large cosmological scale, i.e. the IR cut- 
off L = Rehia). This horizon is the boundary of the 
volume a fixed observer may eventually observe. This 
model has been confronted by cosmic observations ex- 
tensively [l3l-[T5|. for recent results, please see and 
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[17|. In the literature, to the best of our knowledge, only 
the geometry information which includes the luminosity 
distance from SN la, the angular diameter distance 
Da from BAO and the full information of CMB from 
WMAP-7yr were used to constrain this model, for ex- 
amples please see dl] and [13. As is well known, to 
discriminate the cosmological models the geometry in- 
formation is not enough due to the degeneracies between 
model parameters. It means that different cosmological 
models can have the same background evolution history. 
However the dynamical evolution would be very different 
even if they have the same background evolution. Which 
is to say the dynamical evolution is important to break 
the possible degeneracy. 

Thanks to the measurement of the cosmic growth rate 
via the redshift-space distortion (RSD) which relates to 
the evolutionary speed of matter density contrast, now 
one can constrain the evolutions of the density con- 
trast S through /(z)(Tg(z), where f{z) = dhiS/dlna is 
the growth rate of matter and <7s{z) is the rms am- 
plitude of the density contrast at the comoving 8h~^ 
Mpc scale. Here h is the normalized Hubble parame- 
ter Hq — lOO/ikm sec~^Mpc~^. Here we should notice 
that the growth rate of structure f{z) has been used 
to constrain the dark energy model and to investigate 
the growth index in the literature, see [l8| for examples. 
However, the observed values of the growth rate fobs = 
are derived from the redshift space distortion parameter 
/3(z) and the linear bias h{z), where a particular fiducial 
ACDM model is used. It means that the current fobs 
data can only be used to test the consistency of ACDM 
model. This is the weak point of using fobs data points. 
Moreover, the measurements of the linear growth rate are 
degenerate with the bias b or clustering amplitude in the 
power spectra. To remove this weakness, Song & Percival 
proposed to use fas{z) which is almost model indepen- 
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dent and provides good test to dark ener gy m odels even 
without the knowledge of the bias or ag [3. Recently, 
the observed values of f{z)as{z) were provided by the 
2dFGRS [131, WiggleZ SDSS LRG '2^, BOSS 13, 
and 6dFGRS [2411. The latest RSD data points were also 
summarized in [25| . For convenience, we show the data 
points used in this paper in Table HI see also Table 1 of 
Ref. [25]. 



fas{z) Survey and Refs 



0.067 0.42±0.06 

0.17 0.51±0.06 

0.22 0.42±0.07 

0.25 0.39±0.05 

0.37 0.43±0.04 

0.41 0.45±0.04 

0.57 0.43±0.03 

0.6 0.43±0.Q4 

0.78 0.38±0.04 



6dFGRS (2012) [24| 
2dFGRS (2004) [20] 
WiggleZ (2011) [21J 
SDSS LRG (2011) [22] 
SDSS LRG (2011) [22] 
WiggleZ (2011) [21J 
BOSS GMASS (2012) [23] 
WiggleZ (2011) [21J 
WiggleZ (2011) [21] 



TABLE I. Data of fas measured from RSD with the survey 
references. See also Table 1 of Ref. 25] . 



So, the main motivation of this paper is to investigate 
the effect of model parameter c to fcr^{z) and to update 
our previous results by including the current observa- 
tional data of RSD as weU as SN la Union2.1, CMB and 
BAO on constraining the HDE model parameter space. 

This paper is structured as follows. In section [TTl we 
give a very brief review of the HDE model where the radi- 
ation is included and the future event horizon is adopted 
as an IR cut-off. The scalar perturbation evolution equa- 
tions for a spatially flat FRW Universe will also be pre- 
sented. In section IIIIl the constraint methodology and 
results will be presented. We give a summary in section 
lYl 



II. BACKGROUND AND PERTURBATION 
EVOLUTION EQUATIONS 

The energy density of the HDE is written as [1] 



Ph 



(2) 



The Friedmann equation for a spatially flat FRW uni- 
verse reads 

= Hi {Vlr^a-^ + Vlwa-^ + f^coo"^) + ^hH^ : (3) 

where il.^ = pi/3M^iH^ are dimensionless energy den- 
sities for radiation, baryon, cold dark matter and HDE 
respectively. Here the scale factor a has been normalized 



to ao = 1 at present. Combining Eq. ([T|) and Eq. ([2]), 
one obtains the differential equation for fi/j [l^ 



-2nh (1 - n,,) 



( E'{x) 
\E{x) 



(4) 
In a 

This equa- 



where ' denotes the derivative with respect to x 
and E{x) — V^^roe^'^^ + fif,oe^^ + f^coe~ 
tion describes the evolution of dimensionless energy den- 
sity of HDE with the initial condition fJ^o = 1 ~ ^ro ^ 
51f,o — ilco- Via the conservation equation of the HDE 
Ph + 3H{pii + pii) = 0, one has the equation of state 
(EoS) of the HDE 

wi - -1 - i^iln^ - _i _ 2\/r4 

3 din a 3 3c ' 

where the definition Wh — Ph/ Ph is used. 

In this paper, the HDE is taken as a perfect fiuid with 
the EoS ([5]) , then in the synchronous gauge the perturba- 
tion equations of density contrast and velocity divergence 
for the HDE are written as 

Sh = -(1 + wh){eh + b- mP^ - wh)SH, (6) 

2 oph 
e, = - 3c^,,,) + ^-^^esh - kW (7) 

1 + Wh 

following the notations of Ma and Bertschinger [26{ , 
where the definition of the adiabatic sound speed 



Ph 
Ph 



Wh - 



Wh 



snii + wh) 



(8) 



is used. When the EoS of a pure barotropic fluid is neg- 
ative, the imaginary adiabatic sound speed can cause in- 
stability of the perturbations. To overcome this problem, 
one can introduce an entropy perturbation and assume 
a positive or null effective speed of sound. Following the 
work of [27j . the non adiabatic stress or entropy pertur- 
bation can be separated out 



PhTh = 6ph - cl ^^Sph, 



(9) 



which is gauge independent. In the rest frame of HDE, 
the entropy perturbation is specified as 



s.eff 



'^l,ad)^h' 



(10) 



where c^ g^-y is the effective speed of sound. Transforming 
into an arbitrary gauge 



=Sh + 3Hil + Wh) 



fc2 



(11) 



gives a gauge-invariant form for the entropy perturba- 
tions. By using the Eqs ^) (ITUl) and pT|). one can recast 
Eqs. and ^ into 
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K 

2' 



i + Wh 
.2 

-J^eff 1,2, 



5h + m{i + wh)-^ 



1 + 



k Sh - k (Th 



(12) 
(13) 



For the HDE, we assume the shear perturbation ah = 
and the adiabatic initial conditions. Actually, the effec- 
tive speed of sound is another freedom to describe 
the micro scale property of HDE in addition to the EoS 
[28j . And, we should take it as another free model pa- 
rameter. The sound speed determines the sound horizon 
of the fluid via the equation 1^ = Cs^eff/H. The fluid can 
be smooth or cluster below or above the sound horizon 
Is respectively. If the sound speed is smaller, the per- 
turbation of the fluid can be detectable on large scale. 
And in turn the clustering fluid can influence the growth 
of density perturbations of matter, large scale structure 
and evolving gravitational potential which generates the 
integrated Sachs- Wolfe (ISW) effects. However, the au- 
thors of j28| have shown that current data can put no 
significant constraints on the value of the sound speed 
when dark energy is purely a recent phenomenon. For 
the HDE considered in this paper, it is related to the fu- 
ture event horizon and would not cluster. So we assume 



the effective speed of sound 



1 in this work. 



III. 



METHODOLOGY AND CONSTRAINT 
RESULTS 



In our previous work [lg|, we have used the SN la 
Union2, BAO and full information of CMB from WMAP- 
7yr to constrain the model parameter space, where the 
effects of model parameter c to the CMB power spectrum 
were also discussed. In Refs. [l3l, we showed that large 
values of c increase the tails of CMB power spectrum at 
large scale, i.e. I < 10, through the integrated Sachs- 
Wolfe (ISW) effect. Here we will focus on its effects to 
the f<Js{z) caused by the different values of c. At first, 
we modify the CAMB package which is the publicly 
available cod^B for calculating the CMB power spectrum 



to include the HDE. We calculate the values of CTg at 
different redshift for the HDE model. We also write a 
subroutine to calculate the growth rate f{z) for the HDE 
model. The growth rate can be obtained by solving the 
following differential equation [2^ 



din 
where 



1 



2^k{a) + - (1 - Weffia)) ^de{a) 



3(a) 
rjfc(a) 

^deia) 

Weffia) 



D{a) 



= {l + z)D{z), 



~ a^mia)' 

_ VldeH"^ 

= a3[l+i«.„(a)]_f^2(a)' 

In a 

c?ln a'w{a'). 



1 



In 



dg 

din a 

9 = 0(,14) 

(15) 
(16) 

(17) 

(18) 







Here D{a) is the amplitude of the growing mode which 
connects to f{a) via the relation / = d\nD / d\x\a. Fi- 
nally, we can obtain the values of /(Tg(z) at different red- 
shift z. To investigate the effects of c to fa^{z), we bor- 
row and fix the relevant cosmological values from our 
previous results obtained in pL6j but take the model pa- 
rameter c varying in a range. The evolution of fcrs{z) 
with respect to the redshift z for different values of c is 
shown in Figure [T] One can read off that the large values 
of c decrease and increase the values of fas{z) at higher 
and lower redshifts respectively from the Figure [T] It 
clues that the /erg data points favor the values of model 
parameter c in a range of [0.69, 0.9]. However, due to the 
sparseness and relative large error bars of the RSD data 
points, the current data sets of fas{z) may not give a 
much tight constraint to the model parameter space. 



To obtain the model parameter space from currently 
available cosmic observations, we use the Markov Chain 
Monte Carlo (MCMC) method which is efficient in the 
case of more parameters case. We modified the pub- 



licly available cosmoMC packag^ 30] to include the 
likelihood coming from the /(Tg(z). We adopted the 7- 
dimensional parameter space 

P = {wfa, uje, Os, r, c, n„ log[10i%]} (19) 



^ http://camb.mfo/. 



^ http:/ /cosmologist.info/cosmomc/. 
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FIG. 1. The fas{z) v.s. the redshift z for different values of model parameter c (the red dashed line is for c = 0.9, the orange 
thick line is for c — 0.696 and the blue dotted line is for c = 0.5 ), where the other relevant cosmological parameters are fixed 
to their mean values obtained in Ref. [l6l ]. Large values of c decrease and increase the values of fas{z) at higher and lower 
redshifts respectively. The black lines with error bars denote the observed data points as listed in Table IH 



the priors for the model parameters are summarized in 
Table nil Furthermore, the hard-coded prior on the comic 
age lOGyr < to < 20Gyr is also imposed. Also, the phys- 
ical baryon density Wb = flbh'^ = 0.022 ± 0.002 from 
big bang nucleosynthesis and new Hubble constant Hq — 
74.2 ± 3.6kms~^Mpc~^ [s^l are adopted. The pivot scale 
of the initial scalar power spectrum kgo = 0.05Mpc~^ is 
used in this paper. 

The luminosity distance cIl from SN la Uinon2.1 [33l |. 
the angular diameter distance Da and CMB power spec- 
tra from WMAP-7yr are used to fix the background evo- 



lutions. For the details, please see Appendix El 

We ran eight chains on the Computational Cluster for 
Cosmos and stopped sampling when the worst e-values 
[the variance(mean)/mean(variance) of 1/2 chains] R — 
1 was of the order 0.01. The global fitting results are 
summarized in the Table ITIl and the Figurc[2] Comparing 
to our previous result c = 0.696t°:o737±o.l32lo.?90 0' 
we find that SN Union2.1 favors large values of model 
parameter c = O.737l:°;[!*^°t°;}^6+0-320 ^j^g^^ ^-j^g ^iST) 

fcrs(z) is included, the values of model parameter c are 
increased to c = 0.750lJ5;|^^^^i;5 i7^t;5J26 which confirms 
the analysis as shown in Figure [1] 
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Pramctcrs Priors 



Mean with errors without /fjg Best fit without /erg Mean with errors with /erg 



Best fit with /fjg 





[0.005,0.1] 


n rino'7rU.UUUlJi(-|-U. llUlH4-t"U . UUluu 

U.UZZ/ _Q QQQ524-0.00100-0.00151 


0.226 


1 n nnnKAOA-f\ nnii'7_i-n nni^^/i 
'■^■^^^"-0.000549-0.00110-0.00153 


0.0226 




[0.01,0.99] 


nil ri"l"0-00446+0. 00888+0. 0135 
^■-'--'-^-0.00440-0.00863-0.0122 


0.111 


nil r\"l"0' 00478 + 0.00883 + 0.0 145 
U.llU_QQQ4gg_QQQgg2- 0.0127 


0.110 


6 


[0.5, 10] 


_|_0.00264+0.005 12 + 0.00781 
1. (J J95_qqq261-0. 00505-0. 00733 


1.0401 


nr>nA +0.00260+0.00532+0.00758 
1.0o94_QQQ27i -0.00530 -0.00791 


1.0392 


T 


[0.01,0.8] 


_|_Q 00674+0.0255 + 0.0447 
0.0896 _QQQ75g_QQ233_Q_Q3g8 


0.0846 


_|_0 00628+0.0250+0 0466 
0. 0888 _QQQ724 -0.0236 -0.0388 


0.0903 


C 


[0.1, 1.5] 


i_n nsTn_i_n i Qfi_i_n T>n 

^- -0.0826-0.148-0.202 


0.713 


U. 'OU_Q Qggg_Q j^-,3_Q 226 


0.733 


ris 


0.5, 1.5 


U.y / Z_Q g^24_0.0243-0.0370 


0.970 


^■'^^^ ^-0.0131-0.0259-0.0403 


0.970 




[2.4,4] 


^0 0343+0.0690+0 108 
'^■0'9'3_ 0,0341 -0.0669- 0.0940 


3.0730 


r, r.,-^-. +0. 0357+0. 0762+0. 1 14 

3.0766_Q Q3gg_Q QgQQ_Q 


3.0795 


o 




n 71 q+0 0183+0. 0346+0. 0510 
' -^''-0.0176-0.0375-0.0592 


0.719 


n 71 7+0 0171 + 0, 0321 + 0, 0461 
^- ' ^< -0,0174-0,0369-0,0587 


0.717 


Age/Gyr 




1^ qni +0 109+0. 220+0. 314 
10.yul_Q j^Qg_Q 216-0.332 


13.886 


iq Qi p+0. 113+0. 240+0. 354 
^"^■^^"^-0.114-0.221-0.366 


13.923 






n 9Q1 +0.0176+0.0375+0.0595 
-0.0183-0. 0346-0. 0508 


0.281 


n noq + 0. 0174+0. 0369+0. 0587 
"J ■'^^'^-0.0171 -0.0320- 0.0460 


0.283 










n 7fiQ + 0477+0. 0910+0, 120 
U. 'DJ_Q Q4Q5_Q Qg2Q_o,io8 


0.766 






in C/7+11S6+2. 409 + 3. 730 
-'-^■'■"^'-1.219-2.273-3.401 


10.302 


in =70 + 1.189 + 2.335 + 3.794 
lU.0/5_i igg_2 442-3,629 


10.770 


Ha 




fi« 787+1 836+3. 847+5. 602 
oa.iai -1,839-3,080-5.757 


68.927 


en <1 .+1,885 + 3.780 + 5.706 
""•^-'^-1.904-3.795-5.216 


68.451 



TABLE II. The mean values with 1,2,3(t errors and the best fit values of the model parameters and derived cosmological 
parameters, where the WMAP 7-year, SN Union2.1, BAO and RSD fas data sets are used. 



To show the effects of RSD data points /(Tg(z) to con- 
strain the model parameters space, the 2D contour for 
model parameter il„ — c is also plotted in Figure |31 From 
this figure, one can read that the region of flm is shrunk 

I 



I 

when fas{z) data points are employed. But in this case, 
the 1,2,3(7 regions of c are enlarged. And the 2D con- 
tour diagram moves little to the top right corner direction 
when fas{z) data points are included. 



With the mean values listed in the Table [TT] for the case 
of SN-)-BAO-l-CMB+RSD, we plotted the evolutions of 
the EoS of HDE with respect to the redshift z in Figure 
m where the shadows denote the 1, 2, 3cr regions from the 
dark to the light respectively. For calculating the la re- 
gion, we consider the propagation of the errors for w{z) 
and marginalize the other irrelevant model parameters 
by the Fisher matrix analysis [13, UHl • If the other irrel- 
evant model parameters are not marginalized, the error 
bars will be underestimated. The errors are calculated 
by using the covariance matrix Cij of the fitting model 
parameters which is an output of cosmoMC. The errors 
for a function / = f{9) in terms of the variables 6 are 



r 



given via the formula [3 SMS 



E 



5/ 

de. 



df 



dl 

de. 



(20) 

where n is the number of the variables. In our case, / 
would be the EoS w{z] 9i) for HDE. And the variables 9i 
are {Vt^h? .VtjT? ,c) for the HDE model. The correspond- 
ing la errors for w{z) are given by 

Wicr{z) ^ w{z)\g^g±an,, (21) 

where 9 are the mean values of the constrained model 
parameters. For a relative large values of c, the HDE be- 
haves like quintessence at present {'Wh\z=a — —0.971 ± 
0.0777 with la error). In 2a regions, it still has 
broad space to behave like phantom even in the future. 
But in 3cr region, it has the possiblity to behave like 
quintessence. So based on this result, we still do not know 
our Universe will be terminated by a cosmic doomsday 
or not in 3cr region. 



IV. SUIVIIVIARY 



Union2.1 and with the addition of RSD data points of 



In this paper, we updated our previous results obtained 
in Ref. with the replacement of SN Union2 by SN 
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FIG. 2. The ID marginalized distribution on individual parameters and 2D contours with 68% C.L., 95% C.L. and 99.7% C.L. 
by using CMB+BAO+SN+RSD data points. 



f(Js{z). We showed the effects of model parameter c to 
fas{z) by fixing the other relevant model parameters and 
found out that RSD fas{z) data points favor larger val- 
ues of c. But due to the sparseness and relative large er- 
ror bars of the RSD data points, the current data sets of 
/crg(z) cannot give a much tight constraint to the model 
parameter c. A global fitting to the HDE model was per- 
formed by combining the full information of CMB from 
WMAP-7yr, BAO, SN Union2.1, with and without RSD 
fas{z) data sets via the MCMC method. The results 
show that RSD data points fcr%{z) can shrink the model 
parameter space SI™ efficiently as shown in Figure [3] but 
cannot constrain the model parameter c very well. When 
the RSD fcrs{z) data points are added, the 2D contour 
diagram moves little to the top right corner direction on 
the 2D Clm — c plane as shown in Figure [3] It means that 
the RSD fcTsiz) data points favor larger values of c and 
flm ■ It confirms our previous analysis as shown in Figure 

m 

To show the evolution of the EoS with errors for HDE 
with respect to the redshift z, we should margnialize the 



other irrelevant model parameters. If not the error bars 
will be under estimated. We marginalized the other ir- 
relevant model parameters by the Fisher matrix analysis. 
And the evolution of the EoS for HDE in Scr region was 
plotted in Figure |4] by adopting the mean values as shown 
in Table ini In this figure one can see that HDE behaves 
like quintessence at present {wh\z=o = —0.971 ± 0.0777 
with la error). In 2a region, it has a wide region to be- 
have like phantom. But in 3a region, it has possiblities 
to behave like quintessence. Then one still cannot con- 
clude whether the future Universe will terminated by a 
cosmic doomsday or not in 3a region. 
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Appendix A: SN la Union2.1, BAO and CMB 



based on the calculation of 



For the SN la, the Uinon2.1 [33| data sets will be used 
in this paper. The distance modulus is defined as 



IJ,th{z) = 51ogio[rfL(z)] +Mo, 



(Al) 



where dL{z) is the Hubble-free luminosity distance 
HodLiz)/c = HodA{z){l + zf/c, with Hq the Hubble 
constant, and = 42.38 — Slogj^g/i through the re- 
normalized quantity h as Hq = lOO/i km s^^Mpc^^. 
Where c?l(z) is defined as 



dL{z) = (l + z)r(z) 
c 



^{z) = 



dz' 



E{z') 



(A2) 
(A3) 



where E'^{z) = H^{z)/Hq. Additionally, the observed 
distance moduli j-iobs{zi) of SN la at Zi are 



y^obsizt) = mobs{zt) - M, 



(A4) 



where M is their absolute magnitudes. 

For the SN la dataset, the best fit values of the pa- 
rameters Ps can be determined by a likelihood analysis, 



{fJ-obsiZt) - IJ.th{P,Zi)Y 



SN 



{51ogioML(P, Zi)] - mobs{z,) -f Af'}' 
2^ —2 



SN 



where M' = ^.q + M \s a, nuisance parameter which in- 
cludes the absolute magnitude and the parameter h. The 
nuisance parameter M' can be marginalized over analyt- 
ically [38,] as 



/-I-CX3 
exp 
-oo 

resulting to 



B 



C 



dM' 



r=A- — +\n[ — \, 



C 



27r 



(A6) 



with 



8 



-0.4 




FIG. 4. The evolution of EoS for HDE with 1, 2, 3cr shadow regions, where the mean values of the relevant model parameter 
are adopted as listed in the Table [U for the case CMB+BAO+SN+RSD. 



SN 



^ = 51 {^^ogio[dLiP,Zi)] - rriobsizi)} ■ Gov,/ • {blog^oldLiP, zj)] - niobsizj)} , 

SN 

B Covy^ • {5 logio[(iL(P, zj)] - m.obs{zj)} , 

i 

SN 



(A7) 



where Cov,„- is the inverse of covariance matrix vifith or 
without systematic errors. One can find the details in 
Ref. [33II and the web siteEl where the covariance matrix 
with or without systematic errors are included. Relation 
(jA5[) has a minimum at the nuisance parameter value 
M' = B/C, which contains information of the values of 
h and M. Therefore, one can extract the values of h and 
M provided the knowledge of one of them. Finally, the 
expression 



xlAP,B/C)^A-{ByC), 



(A8) 



which coincides to Eq. (|A6p up to a constant, is often 
used in the likelihood analysis jH, HI] . Thus in this case 
the results will not be affected by a flat M' distribution. 
It worths noting that the results will be different with or 



http: / /supernova. lbI.gov/Union/. 



I 

without the systematic errors. In this work, all results 
are obtained with systematic errors. 

For BAO data sets, we used the observational results 
df' from SDSS DR7 and A{z) from WiggleZ [2l|. 
The observed values of df"' are gathered in Table [Xj 



z 




survey and reference 


0.20 


0.1905 ± 0.0061 


SDSS [40] 


0.35 


0.1097 ±0.0036 


SDSS [40] 



TABLE III. The " from SDSS DR7 Uo 



where dz = rs{zd)/ Dv{z), rs{zd) is the comoving sound 
horizon at the baryqn drag epoch, Dv{z) = [(1 + 



zfD\cz/H{z)Yl^ 



iryon ar 

Ell, Si- 



Here Da{z) the angular di- 
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ameter distance which is defined as 

r{z) 



Da{z) = 



l + z 



(A9) 



For WiggleZ data points, one calculates acoustic pa- 
rameter A{z) introduced by Eisenstein et al. [4l| 



(A15) 



For the SDSS DR7 data points, the XsDSsiP) is given ^he observed values of A{z) are gathered in Table H] 



as 



SDSS 



(AlO) 



where C ^ is the inverse covariance matrix 

_ / 30124 -17227\ 
^ I -17227 86977 / 



(All) 



To calculate rs(zd), one needs to know the redshift Zd 
at decoupling epoch and its corresponding sound horizon. 
We obtain the baryon drag epoch redshift Zd numerically 
from the following integration [i^ 



The corresponding XwiqqieZ is given as 



iggle 
WiggleZ 



\Wiggle 



2J 



■{A'^{P.z,)~A'"H z,)lAl(S) 



z A {z) survey and reference 
0.44 0.474 ± 0.034 WiggleZ [21] 
0.60 0.442 ± 0.020 WiggleZ [21J 
0.73 0.424 ± 0.021 WiggleZ [21J 

TABLE IV. The A{z) from WiggleZ [U . 



rvo 

T{Vd) = / dri'id 

drj Xe{z)(TT ^ 

az- — 1 

jQ da R 



(A12) 



where R = 3pb/ipj, ax is the Thomson cross-section and 
Xe{z) is the fraction of free electrons. Then the sound 
horizon is 



Tsizd) 



dr]Cs{l + z). 



(A13) 



where Cs — 1/ \/3(l -I- R) is the sound speed. Also, to 
obtain unbiased parameter and error estimates, we use 
the substitution Hj 



dy—, — Trs\,Zd), 



(A14) 



where dz = rs{zd)/ Dv{z), is evaluated for the fiducial 
cosmology of Ref. |40|] , and Zd is redshift of drag epoch 
obtained by using the fitting formula for the fiducial 
cosmology 



where C ^ is the inverse covariance matrix 

1040.3 -807.5 36.8 
C"^ = I -807.5 3720.3 -1551.9 | (A17) 
336.8 -1551.9 2914.9 

Then the total x%ao from BAO is written as 

xlAoiP) - xloSsiP) + Xw^gglez{P)■ (Al8) 

For the /(TsC^), the X/^gl^) is given 

For CMB data set, the temperature power spectrum 
from WMAP 7- year datcS ii are employed. 

Then one has the total likelihood C oc where 
is given as 

X\P) = xIn{P) + xIao{P) + X).,{P) + xImb(P). 

(A20) 

which is used to get the distribution of the model param- 
eter space. 
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